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We present a method of obtaining the entire unitarity triangle from measurements of B — » nK 
decay rates alone. Electroweak penguin amplitudes are included, and are related to tree operators. 
Discrete ambiguities are removed by comparing solutions with independent experimental data. The 
theoretical uncertainty in this method is rather small, in the range 5-10%. 
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Within the standard model (SM), CP violation is due 
to a complex phase in the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix. It has become standard to 
parametrize this phase information using the unitarity 
triangle, whose apex is given by the CKM parameters 
(p, rf) , and in which the interior (CP-violating) angles are 
known as a, (3 and 7 Q. One of the most important tasks 
in high-energy physics is to measure these quantities, and 
to test whether the SM explanation is correct. 

Many methods have been proposed for getting a, 
and 7 (or, equivalently, p and rf). Most require the mea- 
surement of CP-violating asymmetries in hadronic B de- 
cays 2]. The most promising of these involve decays 
which are dominated by a single decay amplitude (e.g. 
B°(t) — ► J/vjK s , 4>K s )- However, a great many decays 
receive both tree and penguin contributions, with dif- 
ferent weak phases, thus spoiling the cleanliness of the 
methods Q. 

A number of years ago, it was shown that an isospin 
analysis of B — > tttt decays allows one to remove the pen- 
guin "pollution" in B®(t) — > ir + ir~ , so that the CP phase 
a can be measured |4|. Subsequently, Nir and Quinn 
(NQ) showed that this technique could also be applied to 
B — > irK decays, giving another way of extracting a ||. 
(Note that this method gives a with several discrete am- 
biguities 0.) Unfortunately, this analysis neglects elec- 
troweak penguin operators (EWP's), and such operators 
are very important in B — > irK decays 7]. When one 
includes EWP's, the NQ B — > irK analysis fails — one 
cannot obtain weak phase information. (The validity of 
analyses which rely on SU(3) relations between B — > irK 
and B —* irir decays is also affected by EWP's.) 

Recently, it was shown that, by using Fierz transforma- 
tions and SU(3) symmetry, it is possible to relate EWP's 
to tree operators 0,0. In light of this, the B — > nK 
analysis can be resuscitated and improved. As we will 
show, the entire unitarity triangle can be obtained from 
measurements of the B — » irK decay rates alone! In gen- 



eral, the discrete ambiguities can be removed by com- 
parison with experimental data. The theoretical error in 
this method is rather small, in the range 5-10%. 

Using isospin, the B — > ttK amplitudes satisfy a 
quadrilateral relation: 



A+° + V2A 0+ = V2A 00 + A- 



(1) 



where we have defined^' = A(B -> n l K^). The CP- 
conjugate amplitudes A^ satisfy a similar relation (note: 
A +0 corresponds to B~ — » ir~K a , etc.). It is possible to 
express all amplitudes in terms of a number of distinct 
operators. This is equivalent to a description in terms of 
diagrams || . Neglecting the exchange- and annihilation- 
type diagrams, which are expected to be small for dy- 
namical reasons, but including EWP's, there are five di- 
agrams which contribute to B — > irK decays 0, : (1) 
a color-favored tree amplitude T, (2) a color-suppressed 
tree amplitude C, (3) a gluonic penguin amplitude P, 
(4) a color-favored electroweak penguin amplitude P BW , 
and (5) a color-suppressed electroweak penguin ampli- 
tude P ° Jlal . The B — > ttK amplitudes can then be 
written ll 




A+° 
V2A°+ 
V2A 0Q 

A-+ 



P - -P c 

1 q EW 3 



= -P- TeP - Ce l "< - P KW - -Pi 



P — C,p l ~t - P P c 



= -P - Te 11 - -Pi 



(2) 



Here we have explicitly written the weak phase (7), while 
P, T, etc. implicitly include strong phases. To obtain 
the amplitudes A 11 ' for the CP-conjugate processes, one 
simply changes the sign of the weak phases. We have 
assumed that the b — > s penguin contribution P is dom- 
inated by the internal i-quark, so that it has no weak 
phase in the Wolfenstein parametrization 0] (the EWP's 



2 



are known to be dominated by the internal t-quark). In 
this case the amplitudes A +0 and A +0 are identical. 

In the above, the (complex) B — > irK amplitudes are 
written in terms of the six complex theoretical quanti- 
ties P, T, C, P E w, Pew an( i e%1 ■ First, suppose that 
EWP's are absent. In this case it is possible to invert 
the expressions for the amplitudes in order to write the 
theoretical quantities in terms of the magnitudes and rel- 
ative phases of four of the A 13 and A lJ . Now, we can get 
the magnitudes of A 13 and A lJ from measurements of the 
B — > nK branching ratios. However, in order to obtain 
the relative phases, we must fix the ^-quadrilateral and 
the A-quadrilateral, and know their relative orientations. 
In order to do this, we need two additional (real) rela- 
tions involving the A %3 and A 13 . In the absence of EWP's, 
such relations exist. They are: 



A-+ + V2A 00 
V2A 00 + V2A 0+ 



A-+ + V2A 00 , 
V2A 00 + V2A ^ 



(3) 



where A 13 = e 1 A v . (A third relation, not indepen- 
dent, is A +0 + A~+ = A +0 + A~ + .) The first equation 
above indicates that the A- and A-quadrilaterals share a 
common diagonal, the isospin-3/2 amplitude: 



^3/2 = A- 



+ V2A 00 = -(T + cy 



(4) 



Obviously, since the diagonals are common to both 
quadrilaterals, they have the same length. The second re- 
lation is used to determine this length. With this knowl- 
edge, we can fix the A- and the A-quadrilaterals, and de- 
termine their relative orientations. Thus, in the absence 
of EWP's, it is possible to solve for the six theoretical 
quantities: \P\, \T\, \C\, two relative strong phases, and 
7. This is the NQ method 0. 

Unfortunately, in the presence of EWP's, it is no longer 
possible to do this 0, Il2| . In this case there are six 
theoretical quantities, but only five independent ampli- 
tudes in Eq. J5J. Thus, it is impossible to express the 
theoretical quantities in terms of the A 13 and A 13 . (It 
is also straightforward to verify that Eqs. J3J above no 
longer hold if P EW and P£ w are nonzero.) It therefore ap- 
pears impossible to obtain weak phase information from 
B — > irK decays. 

Fortunately, to a good approximation, the EWP's are 
not independent quantities — P EW and P° w can be re- 
lated to T and C. Briefly, the argument goes as follows 
[HI EH The SM effective weak hamiltonian for B — > ttK 
decays is: 



Gi 
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H\ fS = ^= [KfcKUeid + C2O2) - J2 v tbVt* s ^O t ] + h.c. 

v 2 i=3 

(5) 

In the above, 0\ and O2 are (V — A) x (V — A) tree 
operators, while O7-O10 describe the electroweak pen- 
guin operators. O7 and Og have the Lorentz structure 



(V—A) x (V+A), while O g and O w are (V—A) x (V—A). 
However, the Wilson coefficients C7 and cs, which multi- 
ply O7 and 0%, are much smaller than eg and ciq |l4| : 



c 7 = 3.49 x 1(T 4 
c 9 = -9.92 x 1(T 3 . 



c 8 = 3.72 x 1CT 4 
cio = 2.54 x 1(T 3 



(6) 



Thus, the EWP's are approximately given purely by Og 
and Oio- Furthermore, these operators can be Fierz- 
transformed into 0\ and O2, since all have a (V — A) x 
(V — A) structure. Therefore the EWP's are related to 
the tree operators. 

There are two independent relations between EWP's 
and tree operators. Ignoring exchange- and annihilation- 
type diagrams once again, they are given by flU 



P Elv {B A 



~K°) + V2P E 



V (B+ -> 7T°K + ) = 
C 9 + do (T + C) 



P l 



W {B° 

1 Cg - 



d - 

-» TT-K+) + P EW (B' 

dO (T - C) 1 Cg - 



C2 K b v us \ ' 

-» 7T + K°) = 

- do (T + C) 



2 d - c 2 



v: h v u . 



2 d+d \V* b V r 



(7) 



Using the expressions for the B — > ttK amplitudes given 
in Eq. (|2J), these give 

3 c 9 + do 



P, 



P' 



4 ci + c 2 

3 Cg + do 

4 ci 



C2 



R(T- 



do 



C) 



4 ci — c 2 

3 Cg - Cio 



4 ci 



C2 



R(T 
R(T 



C) , 
C) (8) 



where R = \V t * b V ts /V* b V us \. These provide the relations 
between the diagrams P EW and P^ w and the tree opera- 
tors T and C. 

In Refs. [Iollll| . a numerical value is taken for R. But 
in our method, this is not necessary. Instead, we keep 
the general expression 



V t tVts 



V* V 

V,, h Vus 



1 



A 2 VP 5 



(9) 



As we will see below, this allows us to improve consider- 
ably upon the original NQ method. 

The B — > ttK amplitudes are now written in terms of 
the five (complex) theoretical quantities P, T, C, e n and 
\J p 1 + rj 2 . One can invert these expressions to write the 
theoretical parameters in terms of five independent A 13 
and A lJ amplitudes ^5|- However, as discussed earlier, 
in order to determine these parameters, one needs two 
additional (real) relations to fix the two quadrilaterals 
and their relative orientation. There are several ways 
to obtain these. One is to note that, at this stage, e n 
and \J (p- + rf are simply arbitrary complex quantities, 
and are expressed in terms of the A %] and A 13 . However, 
there are physical constraints on these parameters. They 
are: 



0*7 



1 



Im 



2 - 







(10) 
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These provide the relations necessary to fix the relative 
orientations of the two quadrilaterals. 

Note that, in light of the relations between EWP's and 
tree operators, the A 3 / 2 amplitude in Eq. is given by 



.4 



3/2 = -(T + C) 



,11 



3 c 9 + cio 

2 a 



R 



C2 



(11) 



This is one of the diagonals of the ^-quadrilateral. The 
corresponding diagonal in the ^4-quadrilateral is given by 
the above expression, but with 7 — * —7. The relations 
in Eq. ifTUI) imply that |^4 3 / 2 | = \A 3 / 2 | , as was true for 
the case without EWP's. The magnitudes and relative 
phases of the B — > nK amplitudes are therefore obtained 
by measurements of the branching ratios and the con- 
struction of the A- and A-quadrilaterals. This allows us 
to obtain all the theoretical parameters. 

The key point is that there is enough information in the 
B — > irK system to extract the values of seven theoretical 
parameters: the magnitudes of P, T and C, two relative 
strong phases, and two pieces of weak-phase information 
(which we take to be 7 and \J ' p 2 + v/ 2 ) 16]. Note that the 
knowledge of 7 and a/ p 2 + r\ 2 is sufficient to pin down the 
shape of the unitarity triangle. Thus, one can obtain the 
full unitarity triangle (up to discrete ambiguities) from 
measurements of the B — > irK rates alone. 

We now demonstrate numerically how the method 
works. Ideally, we would use current experimental data 
on B — > nK rates. Unfortunately, although the vari- 
ous branching ratios have been measured, no significant 
partial-rate asymmetries have yet been observed , and 
our method requires at least one measurement of direct 
CP violation. We therefore generate values for the "ex- 
perimental measurements" by assuming input values for 
P, T, etc. We choose 

|P| = 1.0 , (5 P = -18.0°, |T|=0.3, 5 T = 2.0° 
\C\ = 0.05 , S = 102.0° , p = 0.18 , t] = 0.38 . (12) 

With these inputs, we find 

\ A +o\ = 1^4+01 = L00 



\A°+\ 


= 0.86 , 


A°+\ 


= 1.00 , 




\A 00 \ 


= 0.62 , 


A 00 \ 


= 0.57 , 




\A~+\ 


= 1.07 , 


A-+ 


= 1.22 . 


(13) 



Here we have taken the values for eg and cio give n in 
Eq. ©, along with c x = 1.144 and c 2 = -0.308 Q. 

Given this "experimental data," we can solve the 
system for our seven theoretical unknowns. Since the 
equations are nonlinear, there will be many discretely- 
ambiguous solutions. For the "data" in Eq. (|13fl . we find 
16 solutions. Half of these yield unitarity triangles which 
point down, i.e. 77 < 0. However, we know from the 
kaon system that rj > 17). We therefore exclude solu- 
tions with 7/ < 0. The 8 remaining solutions are shown 



in Table Q] (Note: we have solved the system for many 
different inputs [Eq. I|12|) ]. In all cases, we find either 
16 or 8 solutions, half of which can be rejected because 
rj < 0.) 

Now, if the SM is correct, there are several constraints 
which these putative solutions must satisfy. First, CP 
violation in B^(t) — > J/ipK s has been measur ed, yield- 



ing a world average of sin 2/3 = 0.736 ± 0.049 18]. Any 
solution in Table [l] which does not give a value for sin 2/3 
in its 3cr range is excluded. Second, the latest 95% c.l. 
range for S = V P 2 + V 2 is 0.356 < S < 0.452 0. An 
acceptable solution must give a value for S in this range. 
Finally, some solutions can be eliminated by making the 
mild theoretical assumption that |P| > \T\ > \C\. (This 
constraint is not essential — we find that solutions which 
do not satisfy this condition generally also violate one of 
the experimental constraints.) 

In the particular case of Tabled the experimental con- 
straints alone eliminate all solutions except (3) (the true 
solution). Indeed, in almost all of the cases we studied, 
in which we varied the inputs in Eq. Ijl2|) . we found that 
only a single solution remained after imposing the con- 
straints. Thus, it is in fact possible to obtain the full 
unitarity triangle from measurements of the B — > irK 
rates alone. 

Note that it is also possible to measure independently 
the indirect CP asymmetry in B®(t) — > ir°K s : 



_ Im (e- 2 ^ 00 M 00 ) 
= |A 00 | \A 00 \ 



(14) 



The knowledge of this quantity will provide a crosscheck 
to the solution(s) found above. If two solutions happen 
to be found, then one can in principle distinguish be- 
tween them through the measurement of A™^" ' . And if 
one finds only a single solution using the above method, 
j^indir famishes an independent check of this solution. 

It is possible that no solution is found which satisfies 
all the constraints and independent measurements. This 
would then be evidence for physics beyond the SM. In- 
deed, the present measurement of the indirect CP asym- 
metry in B a d (t) -> 4>K S may be showing signs of new 
physics: although the BaBar measurement is in agree- 
ment with the SM prediction (within errors), the BELLE 
measurement disagrees at the level of 3.5cr 0]- H this 
discrepancy with the SM is confirmed, this would point 
specifically to new physics in the b — > s penguin ampli- 
tude |20|. Since B — > nK decays also involve b — > s 
penguin diagrams, they would also be affected by this 
new physics. In particular, we would expect to find a 
discrepancy in the values of the parameters of the uni- 
tarity triangle as extracted using the above method and 
in other, independent measurements (e.g. sin 2/3). 

Although there is some theoretical input in this 
method, the uncertainty is rather small. There are three 
sources of theoretical error. First, we ignore annihilation- 
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1^1 


|T| 


\c\ 


(Sp — 5 T ) 


(St — S c ) 


P 


V 


sin 2/3 


s 


a indir 


(1) 


0.96 


0.42 


0.33 


-126.5° 


-28.0° 


-0.56 


0.17 


0.21 


0.59 


-0.41 


(2) 


0.98 


1.97 


1.74 


-21.0° 


-149.2° 


-7.33 


0.99 


0.23 


7.40 


0.31 


(3) 


1.0 


0.3 


0.05 


-20.0° 


-100.0° 


0.18 


0.38 


0.76 


0.42 


-0.80 


(4) 


1.01 


1.90 


0.43 


-31.9° 


-55.0° 


-1.91 


0.18 


0.12 


1.91 


-0.09 


(5) 


1.02 


1.70 


0.55 


-46.7° 


-26.7° 


-0.96 


0.07 


0.07 


0.96 


-0.03 


(6) 


1.02 


1.60 


0.23 


-4.6° 


-8.6° 


-0.88 


0.91 


0.78 


1.27 


-0.48 


(7) 


1.08 


2.05 


0.59 


-5.1° 


-2.4° 


-0.68 


0.37 


0.42 


0.77 


0.24 


(8) 


1.38 


3.12 


1.18 


-5.5° 


-0.6° 


-0.28 


0.06 


0.10 


0.29 


-0.67 



TABLE I: The 8 sets of theoretical parameters which reproduce the "experimental data" of Eq. 1131 . We also give the predicted 
values of each set for sin 2/3, S = y/fP + r? and A™£ ir , the indirect CP asymmetry in B° d (t) -> ir°K s . 



and exch ang e-type diagrams, leading to an error of 
0(1%) BQ. Second, we neglect the Wilson coefficients 
C7 and c 8 compared to eg and C10 , giving an error of about 
4% [see Eq. |J5}]. Finally, one must take SU(3)-breaking 
effects into account in Eq. (JSJ . Ref. 01 estimates such ef- 
fects, and finds them to be roughly 5%. Thus, depending 
on how one adds all the uncertainties, the net theoretical 
error in this method is in the range 5-10%. 

Finally, we must note that experimental errors in the 
B — > ttK branching ratios may make it challenging to im- 
plement this method in practice. However, by performing 
a fit to all experimental data, including constraints from 
sin 2(3 and \J ' p 2 + rj 2 , it should be possible to extract the 
full unitarity triangle. 

In summary, we have presented a method of obtain- 
ing the entire unitarity triangle from measurements of 
B — ► ttK rates alone. It relies on a relation between 
electroweak penguin amplitudes and tree operators. One 
can distinguish among discretely-ambiguous solutions by 
using independent experimental determinations of sin 2/3 
and \J p 2 + vj 2 . The theoretical uncertainty is rather 
small, in the range 5-10%. At present, although B- 
factories have measured the B — > ttK rates, no difference 
between the B and B branching ratios has been observed 
yet. As soon as one observation of direct CP violation in 
B — ► ttK decays is made, it should be possible to extract 
the full unitarity triangle using this method. 

D.L. thanks M. Gronau for helpful conversations, and 
A. Hoecker and L. Roos of the CKMfitter group for the 
latest 95% c.l. range for yj p 2 + rf. This work was finan- 
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